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Mean Coulomb energy of 2d Dyson gas in quadratic potential is examined from combinatorial viewpoint. 
For /3 = 1, we find a recursive relation on mean energy and obtain its exact (finite N) solution in closed form 
in terms of the hypergeometric function 3^2- Using this exact solution, we derive the large- N asymptotic 
expansion of mean energy and show, that this expansion contains half-integer powers of N. 

1 Introduction 

This paper is devoted to the study of a system of N charged particles on the complex plane with complex 
coordinates Zi, which are bound to the origin by attractive Gaussian potential and interact with each other 
by pairwise repulsive Coulomb interactions. In the low energy approximation (neglecting the kinetic energy 
of particles) partition function of such a system can be written as a 2A< r -fold integral 



Z N (fi) = ± f f d 2 Zl ...d 2 z N exp -^Tlz^+^logl^-z/ I (I) 
where (3 determines the strength of Coulomb repulsion. It is conventional to rewrite Z^(/3) as 

Z "(P) = ^n f d 2 z 1 ...d 2 z N H\z l -z J \ 2f3 exp^-Y / H 2 ) (2) 

i<j i 

and to define the correlators (statistical averages in the model) via 

^/(zi,...,zjv)) = ^Z N (P) j d2zi --- d2zN f(zi,---,z N ) Y[\zi - z 3 \ 20 exp (-X)l Zi f) ( 3 ) 



This physical system, known as 2d Dyson gas [1], is quite interesting for a number of reasons. First of all, it 
is expected to undergo a phase transition at certain large value of the Coulomb charge j3 [2, 3]. Second, its 
thermodynamic limit (i.e, the limit of large N) is known to be related to growth processes, such as Laplacian 
growth [4]. The last but not the least, 2d Dyson gas is a typical representative example of /3-ensembles 
[5] , which today attract increasing attention because of their recently established relation [6] to conformal 
field theory and the AGT conjecture [7]. In this context /3-ensembles are used to describe conformal blocks 
of Virasoro and W-algebras and equivalent Nekrasov functions. The parameter (3 is directly related to 
the central charge c in conformal field theory and to Nekrasov parameters e!,e 2 in supcrsymmetric gauge 
theory. All this indicates, that dependence of Zn on j3 is non-trivial, interesting and worth investigating. 
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One of the simplest ways to study this dependence is to consider Zjy((3) as a Taylor series near (3 = 1 



log Z N (f3) = log Z N (l) + ((3- l)E N (l) + (/? 2! 1)2 H N (1) + ... (4) 



and calculate several first coefficients in this expansion as correlators 



E N (l) = (^ l °s\zi-z j \ 2 \ 

\i<3 I 



(5) 



H n( 1 ) = \^^^g\ Zl - z \ 2 \og\z k - Zl \ 2 \ - r^\og\ Zl - z \ 2 \ /^log|z fe -zi\ 2 \ (6) 

\i<j k<l I \i<j I \k<l I 

evaluated at the point (3 = 1. The point /3 = 1 is distinquished, because precisely at this point the system is 
equivalent to the matrix model of complex matrices [8]. Correlators En and Hn have a physical meaning of 
average Coulomb energy and heat, respectively. In this paper, we find the simplest of these coefficients E N 
by combinatorial methods. We find that E N satisfies certain closed recursive relation with respect to N, 
solve this relation explicitly and obtain an exact formula for En- We then derive its large- iV asymptotics, 
which may be of interest in particular physical applications. 

2 Partition function 

Let us note, that a wide class of correlators in 2d model can be explicitly evaluated using a simple rule 

J d 2 z z l z 3 e~ zi = tt % T(I + i) (7) 
One of the simplest examples of such evaluations is the partition function itself at /? = 1: 

^v(l) = ^ J d 2 z 1 ...d 2 z N \{{z l -z J )J\{z l -z J )e^^-Y,\4 2 ) ( 8 ) 

i<j i<j i 

Opening the brackets in the Vandcrmonde determinant and using the rule (7), it is easy to find 



JV JV+l 

^v(i)=n r ( i+fc )' log mi) = E io s r ( i+fc ) ( 9 ) 

fc=0 fc=l 



From the point of view of the Taylor expansion (4), this provides an exact formula for the Oth term. The 
rest of this paper is devoted to a similar (though technically more involved) calculation of the 1st term of 
the expansion (4). Higher order terms would be also very interesting to compute; however, they are more 
complicated and remain beyond the scope of present paper. They will be considered elsewhere. 
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3 Mean energy 

We find it convenient to evaluate the mean energy 



En{1) = n N Z N (l) J d2z ^-- d2zN [ E lQ gl^ ~ z i\ 2 1 IlN-^| 2ex P(-El^l ) ( 10 ) 



using the replica limit 



E N {1) 



lim 

e-s-0 



eAf(e) 



,2e 



(11) 



where the correlator 



e N {e)^- w ^j^(d 2 z l ...d 2 z N y^\z t -z \ 2t \ ]J \ Zi - Zj \ 2 cxp ( - ^ |z;| 2 ) (12) 

can be easily calculated for natural e with the help of (7). We present this calculation here in some detail. 
Since the integrand is symmetric under all permutations of zi,...,Zn, we have 

e w( £ ) = m / d2z i--- rf2z w|zi-2: 2 | 2£ ]J|z J -z 3 | 2 cxp(-^|z l | 2 ) (13) 

' i<j i 



Expanding the Vandermonde determinant, we find 



m 



-N z l 



=3n 



(14) 



where the two sums over (ii, . . . , ijv) and (ji , • • • , jjv) are not written explicitly but assumed. Both sums 
are taken over permutations of (0, . . . , N — 1), and £ is the sign of permutation. Using (14) we obtain 

ew(e) = 2t!^Z n (1) J d ^ Zl ■ ' ■ ^ ZN ' Zl ~ Z2 ^ z ^ ' ' ' z n z i ' ' ' z n cxp ( ~ zZ ( 15 ) 



Integrating out variables zs,...,Zn with the help of integration rule (7) we find 



ejv(e) 



N(N-l) 



N 

II' ^ / 

k=3 J 



d 2 Zl d 2 z 2 \ Zl - z 2 \ 2e z^z^zfzi 2 exp 



\Zl\ 



R2 



(16) 
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With the help of the same rule (7), the partition function Zn(1) can be calculated: 
Z N {1) = I d2z i ■ ■ -d 2z N \\ \zi - Zj^cxp ( - ^2 \ z i\ 2 ) = 

= ;V / d 2 Zl . . . d 2 z N z{- . . . z^z? . . . ztf exp ( - ^ |z 4 f) = [] fc! 

Substituting (17) into (16), we obtain 

AT (AT — 1) £i 1 i 2 i 3 ...i N £j 1 j 2 i3---iN^- ■ ■ ■ ^N* / ,2 ,2, . „ ,2e n i 2 -j! -j 2 / I |2 | |2\ 

ejv(ej — ^-"2 — — / a z Y a z 2 \z\ — z 2 \ z 1 z 2 z 1 z 2 exp ^ - \z\\ - \z 2 \ j 

" (17) 

After minor simplification, this expression takes form 

e ^( £ ) 4 E \\ I d 2 zid 2 z 2 \zi - z 2 \ 2e (z a x z\z a x z\ - zlz^zi) exp ( - \ Zl \ 2 - |z 2 f) (18) 

n 0<a<b<N °" ' •* 

At natural e, this correlator can be calculated by opening the brackets 



i* - ^ - E o i!(e l , jit— ), (-D i+ ^r4r^ a d9) 

and using the same integration rule (7). In this way we find 

n _ V /y e! 2 ^° (-i)q-fc e! e! \ (e - » + «)!(z + 6)! 

I Z !2( e _,)i2 i!( c -i)! (i + 6_ a )!(e-i-6 + a)! i a!6! 1 j 

Note, that for each particular natural A" and e formula (20) is just a finite sum, which can be evaluated 
exactly. Given this combinatorial formula, one straightforwardly proves that ejv(e) satisfies a recursion 



(N + l)(e + 2A^ + 3)(e + 2A^ + 4)e N (e) - 

- 2(12 + SAN + 26N 2 + 6N 3 + 5e + l(We + AN 2 e + e 2 + Ne 2 )e N+1 {e) + 
+ (N + 1)(12N 2 + 38N + ANe + e 2 + 5e + 24)e N+2 {e) - 

-4(N + 2)(N + l) 2 e N+3 (e) = (21) 





e — 1 


e — z 


e — o 


€ — 4 


e — o 


N = 2 


4 


24 


192 


1920 


23040 


N = 3 


15 


114 


1152 


14400 


213120 


N = 4 


36 


332 


4056 


60720 


1064640 


N = 5 


70 


760 


10890 


189720 


3838200 


AT = 6 


120 


1500 


24660 


489744 


11218320 


N = 7 


189 


2674 


49602 


1105608 


28268520 



Figure 1: Several first values of e N (e), calculated with the help of (20). 

The proof is just by substitution. This is very useful: expanding now eAr(e) in powers of e 

e N (e) = N{N ~ 1) +eE N (l) + ... (22) 
we find a recursion, satisfied by the mean Coulomb energy: 

2(2iV + 3)(JV + 2)(N + l)E N (l) - 

- 4(JV + 2)(W 2 + 7N + 3)E N+1 (1) + 
+ 2(N + 1)(6N 2 + 19N + 12)E N+2 {1) - 

- A(N + 2)(N + l) 2 E N+3 (l) + (3JV + 5)(JV + 1) = (23) 
The general 3-parametric solution of this linear 3-rd order difference equation is 

where tf(AT) = EaLi V fc - E q- ( 24 ) 

can be easily proven by direct substitution. From the initial conditions 
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f E 2 (l) = l- 7 + ln2 



£ 3 (1) 



29 



37 + 31n2 



793 

£ 4 (l) = 7^-67 + 61n2 



we find the following values of the constants C\ , c 2 , c 3 : 



13 1 1, „ 

a = — h -7 — In 2 

1 8 2' 2 



9 1 , „ 
C2 = -8 + 2 ln2 



, c 3 = 1 



Therefore, we obtain an exact formula for mean Coulomb energy at finite N: 



E N (1) = ±N 2 *(N) + ±+C^ + ^-^\n2\N+(~ + ^\n2) Y f 



^ NjN-l) 1 T { k+ 2 
^ fc(fc-l) (fc + 1)! /3 



The remaining sum can be evaluated in terms of the hypergeometric function 3 F 2 : 



^ fc(fc-l) (fc + 1)! r /3 



N(N-l) (1 



In 2 - 



r(iv +2 



r (n + 2) r 



if 2 



1,^-1,^ + 3/2 

7V + 2,7V + 1 



Using this formula, we finally obtain a relatively simple answer: 



N^(N) JV 2 3JV 1 JV 7 r (TV + 3/2) 
^iv(l) - 2 T + X + 4 + ^~ r(7V + 2)r(3/2) 3 ^ 2 



1,7V- l.JV + 3/2 
7V + 2,7V + 1 



(25) 



If one rescales partition function 



Zn(P) 



/^zi^.^ZTYj^l^-Zj^exp^-TV^Izjl^ = ( 



, s 0N(N-1)+2N 

im) ZM (26) 



Figure 2: Mean Coulomb energy (25) for small N. 



then one also needs to shift the energy: 

E N (l)=E N (l)-N(N-l)/2hiN (27) 



The values of -Ejv(l), calculated with this formula, coincide with Monte-Carlo results, obtained in [9]: 



N 


Exact answer via (27) 


Its evaluation to 20 digits 


Monte-Carlo result 


Monte-Carlo error 




5 


5831 10 7 + 10 In ( 
384 ' \ 


I) 




0.249831865576454075 


0.25013 


0.00030 




10 


(...) -45 7 + 45 In ^ 


h) 




-7.766078379434240251 


-7.76667 


0.00065 




25 


(...)- 300 7 + 300 In ( 


' 2 A 




-99.84816319738693935 


-99.85156 


0.00329 




50 


(...)- 12257 + 1225 In 


{- 
\50 


) 


-492.8773199345421588 


-492.88918 


0.00702 




75 


(...)- 27757 + 2775 In 


(1 


) 


-1191.5028305194359964 


-1191.49356 


0.00974 

















where (. . .) stand for certain rational numbers, which we do not write explicitly due to space limitations. 
This good agreement with numeric calculations supports validity of the main formula (25). 
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4 Large N expansion of the mean energy 

Let us now obtain the large- N expansion of the mean energy (25) . Since 



T k 



fc(fc-l) (fc + 1)! 



^/7T 4V7T 64V7T 



we have 



V 1 1 ^ = const - — N~ 3 / 2 + 7 TV" 5 / 2 - 391 N~ 7 ' 2 + 

2^k(k-l)(k + l)\ T (3\ C ° nbt 3^ + 10^ 672^ 



where the constant is easy to find: 



C ° nSt = E fc(fc _ 1)(fc + 1)I r ^ = 8 " 2 ^ 



Consequently, 



N-i „,„ , r | fc+ 3 



£^ fc(fc-l) (fc + l)! /3\ v ' \8 2 J 37^ 307^ 33607^ 



k=2 v y \ 1 /• p 
Similarly 

N 2 ,„s „ AT 1 1 „ , 



*(iV) = -AT 2 In AT + A^ 2 A^ 4 + ... (28) 

2 v 7 2 4 24 240 504 v ' 



After all cancelations and simplifications, we find the large- Af expansion of (25) 



^ 1 r2l Ar iV 2 (I 1 \ r 5 AN 1 / 2 N-V* 107N- 3 / 2 N~ 2 , on , 

-Eiv(l) = — A^ lnAf hi — I — 7 1 AH p- =^ H h . . . (29) 

ww 2 4 V 2 2 7 24 30F 3O0F 33600? 240 v ' 



which contains, except for the leading logarithmic asymptotics, integer and half-integer powers of AT. 
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5 Conclusion 



In this paper, we have obtained an exact finite ./V formula (25) for the mean energy of j3 — 1 2d Dyson gas 
and its large N expansion (29). This mean energy describes the 2d Dyson gas partition function Zn((3) in 
the first order of the (/3 — 1) expansion. It would be interesting to extend these results to higher orders of 
the (j3 — 1) expansion. In this direction, the next simplest quantity after the mean energy is the mean heat 
Hpf(l). In analogue with E^(l), we expect ffjv(l) to satisfy some recursive relation in N. This relation, 
as well as its exact solution, remains to be found and investigated. It is also desirable to find a conceptual 
explanation for appearance of such recursive relations. Probably, the origin of these recursive relations lies 
in the integrable structure [10] of the complex matrix model. 
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